The role of hyperfine coupling in magnetic and quadrupolar ordering of Pr3Pd2oSi6 
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We study the ternary clathrate Pr3Pd2oSi6in specific heat and AC-susceptibihty measurements 
on a high-quality single crystal, distinguishing antiferromagnetic (AFM) and antiferroquadrupolar 
(AFQ) ordering on two sublattices of inequivalent Pr sites. The specific heat shows the direct 
involvement of nuclear spin degrees of freedom in the AFM ordering, which is well supported by 
our calculation of the hyperfine level scheme without adjustable parameters. Pr3Pd2oSi6 is therefore 
one of the rare materials where the nuclear moments are involved in the formation of the magnetic 
ground state. 

PACS numbers: 82.75.-z, 31.30.Gs, 75.30.Kz, 75.30.Cr, 75.30.Sg 
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INTRODUCTION 

The ground state of solids is commonly described with- 
out referring explicitly to the nuclear degrees of freedom. 
In lanthanide compounds there are, however, rare exam- 
ples where nuclear moments directly influence the mag- 
netic ground states. This may happen when the dipole 
hyperfine (HF) coupling A of nuclear spin I and total an- 
gular momentum J of 4/ electrons is sufficiently large. 
In Pr compounds, where A has the exceptionally high 
value of 52 mK [l|, Q , the HF interaction becomes even 
decisive for the magnetic order of the ground state, if the 
crystalline electric field (CEF) effect favors a 4/ singlet. 
For example, in pure Pr metal incommensurate induced 
antiferromagnetism requires the assistance of the large 
HF coupling [1, Ql . HF interactions also play an^essen- 
tial role in the singlet ground state system PrCu2 [5] , and 
they were found to induce magnetic ordering of a S* = 1/2 
quantum magnet on a Kagome lattice Importantly, 
the magnetic transition temperatures of such coupled 
systems are usually on the order of a few milli-Kelvin, 
below the base temperature achievable in most labora- 
tories. Here we provide evidence that in Pr3Pd2oSi6 not 
just the electronic but also the nuclear degrees of freedom 
are involved in a magnetic transition at T/vi — 0.14 K, a 
temperature that is easily reached by a standard dilu- 
tion refrigerator. This material is therefore an excellent 
candidate to investigate the role of HF coupling in the 
formation of magnetic ground states. 

Furthermore, Pr3Pd2oSi6 belongs to a new class of 
ternary clathrates R3Pd2oX6 (R = rare earth, X — Si,Ge) 
that has attracted considerable interest after the discov- 
eries of quadrupolar (QP) ordering and thermally ac- 
tivated rattling Q of Ce ions in Ce3Pd2oGe6. R atoms 
in these cage compounds reside on two sublattices of in- 
equivalent sites: a face centered cubic lattice of 4a sites 
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FIG. 1. (Color online) Cage structure and HF level scheme 
for Pr 4a {Oh symmetry) and 8c sites {Td symmetry). The 
energy levels are calculated at Ji" = 0, using the CEF Hamil- 
tonian obtained from INS [l^. Inset graphs: calculated H- 
dependence of the lowest two CEF levels for iif || [110]. The 
line thickness indicates the scale of HF splitting. 



with cages of cubic Oh symmetry (Fig. (TJ left), and a 
simple cubic lattice of 8c sites, where the cages have cu- 
bic Td symmetry (Fig. [TJ right). The occupancy ratio of 
8c:4a is 2:1. Depending on the 8c and 4a CEF ground 
states, a rich variety of multipolar orders may appear 
at low temperatures, as witnessed by measurements of 
the elastic constants [10Ml2|. as well as thermodynamic 
and transport studies [9|. The Pr3Pd2oSi6 clathrate in- 
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FIG. 2. (Color online) (a); Specific heat C{T) of 
Pr3Pd2oSi6 at B = T, measured for SC (black symbols) and 
PC samples (gray symbols). While the SC measurements 
clearly resolve two peaks at Tq = 0.21 K and Tni ~ 0.14 
K, these features are absent in the PC data. Instead, a broad 
peak is observed around 0.6 K. Solid lines show the calculated 
specific heat for the 4a (red) and 8c sites (blue) , and the sum 
of the two (green) . Non-magnetic contributions (gray line) are 
estimated from La3Pd2oSi6 [l^. (b): Calculated entropy S 
(green line), compared to the entropy obtained from the mea- 
sured C/T (black). The additional entropy release of 2i?ln2 
between 0.07 and 2 K indicates fully developed ordering of 
the 8c sites, lifting the two-fold degeneracy of the Fa ground 
state. 



vestigated here is particularly interesting, since the CEF 
ground state entails the possibility of both magnetic and 
QP order. Under the cubic symmetry of the CEF Hamil- 
tonian TLcef, the 9- fold Pr 4/ multiple! (J = 4) splits 
into two triplets (F4 and F5), one F3 doublet and one 
Fi singlet. Inelastic neutron scattering (INS) has 
shown different CEF ground states for the two Pr sites: 
a Fs triplet on the 4a sites and a nonmagnetic F3 doublet 
on the 8c sites, with electric quadrupoles as the lowest- 
order moments. In this Letter, we study both the mag- 
netic and quadrupolar ordering of the two Pr sublattices 
in Pr3Pd2oSi6, and we show how the large on-site HE 
coupling between the 4/ and the i^^Pr nuclear moments 
influences these ordering transitions. 



SINGLE-ION CALCULATIONS 

We calculate the HE level scheme, as described in 
Ref. [3, by adding the on-site HE dipole interaction AI J 
between the 4/ total angular momentum and the i**ipj. 
nuclear spin (/ == 5/2) to the experimental CEF Hamil- 
tonian from Ref. 13: 

U = HcEF + AI ■ J + [gj^isJ - gN^J'NI) ■ Mo-ff- (1) 

An external magnetic field ^qH induces electronic and 
nuclear Zeeman splitting (gj^BJ — gN^J■NI) ■ ^J'oH , where 
gj — 4/5, gN = 1.72, and /xb or fi^ represents the Bohr 
or nuclear magneton, respectively. Using the coupling 
constant A — 52 mK from Refs. [H we calculate the 
eigenenergies of the Hamiltonian ([T]) shown in Fig. [T] For 
the 4a site, where the CEF ground state is a F5 triplet 
with a nonzero magnetic dipole moment, the large on- 
site HE dipole interaction significantly affects the zero 
field energy spectrum and causes a reconstruction of the 
F5 triplet into a quartet at K, a sextet at 0.323 K, 
and an octet at 0.778 K, with additional splitting on the 
order of mK. The non- magnetic F3 doublet on the 8c-site 
shows only small HE splitting of 4 - 5 mK due to second 
order effects in A. The magnetic field dependence of 
the lowest energy levels is dominated by the electronic 
Zeeman splitting of the F5 triplets. Large fields applied 
along the [110]-direction lead to a splitting of the non- 
Kramers F3 doublets above ^qH « 3 T for the 4a sites 
and w 6 T for the 8c sites. 



EXPERIMENTAL RESULTS 

To distinguish magnetic and non-magnetic phases in 
Pr3Pd2oSi6, and to obtain a complete H — T-phase dia- 
gram at low temperatures, we combined measurements of 
the specific heat C and magnetic AC-susceptibility xac- 
All measurements were made on the same Pr3Pd2oSi6 
single crystal (SC), which was fabricated by floating zone 
melting using an image furnace with two ellipsoidal mir- 
rors. The specific heat was measured with a semia- 
diabatic compensated heat-pulse method [l5|, using a 
'^He/'^He dilution refrigerator between 0.07 and 4K, and 
a Quantum Design PPMS equipped with a '^He-insert 
between 0.3 and 10 K. The xac measurements between 
0.02 and 4 K were carried out in a Kelvinox 400 •^He/'* He- 
dilution refrigerator, with a modulation field of 15 f/T 
at 113.7Hz. Absolute values of xac were determined 
by matching the data above 2 K to calibration measure- 
ments in a Quantum Design SQUID VSM. 

Fig. [2] shows the specific heat C measured at zero 
magnetic field. We observe two distinct peaks at Tq = 
0.21 K and Tjvi — 0.14 K, consistent with phase tran- 
sitions found iri previous ultrasound measurements of 
Pr3Pd2oSi6 SC [lllllll- These transitions were not ob- 
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FIG. 3. (Color online) (a) Calculations of Cpr = Cp, (4a) + 
Cpr{S:c) for II [110]. (b)T-dependent specific heat C, mea- 
sured at various fields OT < ^oH < 9T. For clarity, the 
curves in panels (a) and (b) are shifted by +4 J/Kmol per T 
relative to H = 0. Tjvi and Tq indicate the critical temper- 
atures for phase transitions at low fields. Comparing panels 
(a) and (b), further peaks in C{T) and x' {T) at higher fields 
are identified as Schottky peaks from Zeeman splitting of the 
electronic and HF ground state, (c) AC-susceptibility x {T) 
at dc-fields between and 1.2 T. The peak at Tmi appears in 
x'(r) while the peak at Tq is absent. An additional transition 
is found at Tn2 = 0.06 K. 
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FIG. 4. (Color online) Real part x! {T) and imaginary 
part x!'{T) of the magnetic AC-susceptibility, at various 
dc fields < ^loH < 105 mT. x'{T) curves are offset by 
-1-0.4 X lO'^^m^/mol per 10 mT relative to H = 0. Green 
and blue arrows indicate the transitions at Tjvi and Tn2, re- 
spectively. Inset: four-point resistance measurement of the 
same crystal at OT, showing a change of slope around Tn2- 



served in specific heat data of polycrystalline (PC) sam- 
ples, prepared by arc-melting and annealing [l6| . The 
3-20-6 phase is known to possess a very narrow ho- 
mogeneity range [l3|, causing a strong sample depen- 
dence of physical properties. Our specific heat data show 
that high-quality SC are necessary to study the low- 
temperature phases in Pr3Pd2oSi6. To quantify the spe- 
cific heat related to the ordering of both Pr sites, all non- 
magnetic contributions must be subtracted. These are 
estimated from reference measurements of the isostruc- 
tural non-magnetic compound La3Pd2oSi6 [l6| (gray line 
in Fig. [2]), and turn out to be negligible in the region of 
interest T < 2 K. Next, we compare the measured specific 
heat to the single- ion specific heat Cpr (Fig. [21 green line) 
calculated for the two Pr sites using the HF level schemes 
in Fig. [1] We find that the low-temperature specific heat 
is greatly enhanced by excitation of nuclear degrees of 
freedom. Integrating the measured C/T over T (Fig. [51 
bottom panel), we obtain AS' = 24.05 J/Kmol for the 
entropy between 0.07 and 2K. The degeneracy of the Fs 
and Fs CEF ground states at the 4a and 8c sites can only 
account for In 3 -I- 2i? In 2 = 20.66 J/Kmol, if both sites 
order completely. Without the involvement of nuclear 
degrees of freedom, the measured AS would therefore be 
in st rong contrast to the CEF ground state obtained from 
INS [13[. However, our calculation shown in Fig. [2] pre- 
dicts a large specific heat peak due to the thermalization 
of HF multiplets on the 4a sites, which coincides with the 
low-T phase transitions. The expected entropy release 
from the thermal occupation of the 4a site HF sextet at 
0.323 K and the octet at 0.778 K is i?(ln(18) - ln(4)) = 
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12.51 J/Kmol. Subtracting this value from the measured 
AS* leaves an entropy of 11.54 J/Kmol, which is very 
close to i?ln(4) = 2i?ln(2) = 11.53 J/Kmol. 

The additional entropy release could either be at- 
tributed to ordering of the HF ground state quartet of 
the 4a sites {AS = i?ln(4)), or to fully developed order- 
ing of the doubly occupied 8c sites, lifting the two-fold 
degeneracy of the Fa ground state(A5' = 2i?ln(2)). XAC- 
measurements can distinguish these magnetic or non- 
magnetic types of ordering. Fig.[3]compares T-dependent 
specific heat C (T) (panel b) and the real part x' of Xac 
(panel c) to calculations of Schottky anomalies in C{T) 
(panel a), based on the previously calculated single- ion 
level scheme (cf. Fig. [T]). The peak positions in C{T) 
agree well with the calculations. Only the field-induced 
splitting of the F3 doublets is significantly underesti- 
mated by the numerical results, where the measurement 
at 7T closely resembles the calculated curve at 10 T. We 
observe a common peak in C and x' &t T^i = 0.14 K, 
whereas the second peak at Tq = 0.21 K only appears 
in the specific heat for all the applied magnetic fields. 
This implies ordering of the non-magnetic 8c sites at Tq , 
where the lowest-order moments are electric quadrupoles, 
and attributes the entropy release AiS* = 2i?ln(2) to 
fully developed QP order on the 8c sublattice. While 
the F5 ground state of the 4a sites would in principle 
support non-magnetic QP order as well, the transition 
would be revealed by the susceptibility at finite fields, 
due to a field-induced mixing of dipolar and QP mo- 
ments 0, [3, [3| ■ Therefore, the non-magnetic transition 
at Tq is interpreted as QP ordering of the 8c sites, and 
the magnetic transition at Tni as AFM ordering of the 
4a sites. The HF ground state of the 4a sites remains 
effectively four-fold degenerate down to 0.07 K. 

In addition to the transitions at Tq and T/vi, we ob- 
serve a third feature at T]y2 — 0.06 K, presented in more 
detail in Fig. |4l A drop of x' below Tn2 = 0.06 K is ac- 
companied by a peak of the imaginary susceptibility x", 
which indicates dissipative processes. With increasing H, 
T/V2 shifts slightly to lower temperatures and the anoma- 
lies in both x'{T) and x"(T) gradually disappear. A four- 
point resistance measurement (Fig. |4l inset) shows weak 
T-dependence with a slight increase of « 5 fj,fl between 
0.1 and 0.02 K, and a change of slope occurs around TN2. 
Without further measurements it is difficult to establish 
the nature of this transition. 

Fig. [5] shows the phase diagram derived from the T- 
and iJ-dependence of features in C and x' ■ The high- 
est temperature phase boundary agrees very well with 
previous results from measurements of the elastic con- 
stants [l2|- From the shape of the phase boundary and 
symmetry arguments, this transition was interpreted as 
AFQ ordering of O2 = ^i^j"^ - - Jy) moments at 
the 8c sites, where the T3 ground state does not support 
a magnetic dipolar order parameter 12|. Our observa- 
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FIG. 5. (Color online) Phase diagram of Pr3Pd2oSi6 derived 
from specific heat (triangles) and AC-susceptibility XAc(r) 
(circles) and xac{H^ (star symbol). Results from Ref. [121] 
(diamonds) are added for comparison. We identify an AFQ 
ordered phase (magenta) of the Pr 8c sites, followed by AFM 
ordering of the 4a sites (green) below 0.14 K. The additional 
phase below Tn2 ~ 0.06 K (blue) is yet unknown. 



consistent with full ordering of the 8c sublattice confirm 
this interpretation. The phase below T/vi = 0.14K is at- 
tributed to AFM ordering of the 4a sites. The newly dis- 
covered third transition below T]\[2 = 0.06 K must involve 
magnetic moments of the 4a sites as well, since the tran- 
sition shows up in the zero field susceptibility. A first- 
order incommensurate-to-commensurate AFM transition 
is one possible scenario, which would be compatible with 
both the finite x" and the change of slope in the electrical 
resistance, showing only weak T-dependence below Tn2- 
The ordering could also involve the effective magnetic 
moment of the 4a-site ground state quartet (pseudospin 
3/2 [l3|)) which is still degenerate above T7V2 — 0.06 K. 



CONCLUSIONS 

We have shown how specific heat and AC-susceptibility 
measurements can be effectively combined to study mag- 
netic and non-magnetic multipolar phases in Pr3Pd2oSi6. 
We confirmed previously reported AFQ and AFM order- 
ing and discovered an additional magnetic phase tran- 
sition at 0.06 K. Our results further demonstrate that 
Pr3Pd2oSi6 is another one of the rare examples where nu- 
clear degrees of freedom are directly involved in the order- 
ing of electronic 4/ multipoles. Since the calculated HF 
peak in the specific heat of the 4a sites coincides exactly 
with the experimentally observed AFM phase transition, 
it is likely that the HF interaction plays an essential role 
in the AFM order. The present theoretical analysis of 
specific heat and entropy is limited to single ion CEF 
potentials and HF interactions. A more realistic descrip- 
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tion should include inter-site dipolar and quadrupolar in- 
teractions and their corresponding molecular fields po| . 
This will give a deeper insight on how the large Pr HF 
coupling influences the thermodynamics of magnetic and 
quadrupolar order in the Pr clathrate. 
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